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ABSTRACT 

The Heterodyne Instrument for the Far Infrared (HIFI) onboard the Herschel Space Observatory allows the first observations of light diatomic 
molecules at high spectral resolution and in multiple transitions. Here, we report deep integrations using HIFI in different lines of hydrides 
towards the high-mass star forming region AFGL 2591. Detected are CH, CH+, NH, 0H+, H2O*, while NH+ and SH+ have not been detected. All 
molecules except for CH and CH^ are seen in absorption with low excitation temperatures and at velocities different from the systemic velocity 
of the protostellar envelope. Surprisingly, the CH(//r/. = 3/22- - l/2i_+) and CH^(y = 1-0, 7 = 2-1) lines are detected in emission at the 
systemic velocity. We can assign the absorption features to a foreground cloud and an outflow lobe, while the CH and CH^ emission stems from 
the envelope. The observed abundance and excitation of CH and CH* can be explained in the scenario of FUV irradiated outflow walls, where 
a cavity etched out by the outflow allows protostellar FUV photons to irradiate and heat the envelope at larger distances driving the chemical 
reactions that produce these molecules. 
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1. Introduction 

The simplest constituents of interstellar chemistry, light diatomic 
hydrides and their ions (XH and XH+ with X=0, H, C and N), 
have so far not been studied thoroughly. Since their rotational 
ground state lines lie mostly outside atmospheric windows, they 
cannot be observed from ground. However, they are key-species 
of the chemical network. For example, most water at high tem- 
perature (T > 250 K) is formed by the reaction of OH with 
H2. Moreover, once their chemistry and excitation is well under- 
stood, t hey may be valuable tracers of warm and partlv i onized 
gas (e.g. lSternberg & Dalparnoll995l;lBruderer et al.l20l"0l) . This 
is of particular interest in the context of star formation, as the 
radiation of the forming star may heat and ionize the envelope 
giving birth to it. 



Chemistry and excitation of hydrides is characterized by 
high activation energies and high critical densities. For example, 
the only formation route of CH* is via the highly endothermic 
reaction C+ + H2 + 4640 K ^ CH+ + H. The critical densit>Q 
for many hydrides is very high. For example the CH*(7 =1-0) 
transition has a critical density of ~ 5 x 10^ cm~^ for colli- 
sions with H2 at 500 K (iHammami et all 1200^ . This very re- 
active ion may h owever be de stroyed rather than excited in col- 
lisions with H2 (Black 1998). Thus, the excitation of the ions 
studied here may be governed by the initial state in which they 
are formed rather than by equilibrium with collisions or the ra- 
diation field. The situation of considerable abundance only at 
high temperature (T » 100 K), high critical density and quick 
destruction in collisions with H2 is similar to CO*, studied by 
Stauber & Bruderer (2009) . Finally, if a molecule is abundant, 
but not excited, it may still act as an absorber of a background 
continuum source and shows up as an absorption line. 



Send offprint requests to: Simon Bruderer, 
e-mail: simonbr@astro . phys . ethz . ch 

* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



' The density, at which coUisional deexcitation is comparable to 
spontaneous decay. 
^ Scaled for the reduced masses of the H2-CH* system. 
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The Heterodyne Instrument for the Far Infrared (HIFI; 
Ide Graauw et alj2 010^ onboard the Herschel Space Observatory 
allows for the first time to carry out spectrally resolved observa- 
tions of different transitions of hydrides and thus to disentangle 
moving regions along the line of sight and to study fine/hyperfine 
structure components. 

AFGL 2591 is a well studie d high-mass star forming re- 
gion at a distance of ~ 1 kpc (e.g. Ivan derTak et al.i 199 9). The 
sour ce is relatively nearb y and its geometry thus well known 
(e.g. lPreibisch et al]|2003h . There is evidence for a cavity along 
the outflow that allows FUV photons to escape to longer d is- 
tances and irradiate the outflow walls jBruderer et al.ll2009alfa) . 
In these high-density PDRs, the abundances of hydrides can be 
enhanced by several orders of magnitude. Chemical models of 
this scenario predict fluxes for the CH^(y =1-0) line of up to 
16 K km s ' for AFGL 2591, while those from a spherical model 
are fa r below the detection limit of Herschel (iBruderer et alJ 

l2oToh . 



2. Observations and data reduction 

YiWl-Herschel observations of AFGL 2591 

(q'20()0 = 20h29™24^;87, ^jooo = +40°11'19.5") have been 
carried out as part of the Water In Star forming regions with 
Herschel (WISH; van Dishoeck et al., in prep.) guaranteed 
time key program. Target lines are given i n Table [H Molecula r 
parameters are obtai ned from the CD MS (iMiiUer et al.l l2005h. 
except for NH+ dHubers et al.1 l2009h and H9O+ (iMiirtz et alJ 
[1998). Dual beam switch mode with fast chopping was used. 
To disentangle lines from the upper and lower side band, the 
local oscillator frequency was slightly shifted after half of the 
integration time. Data were converted to r,nb using a main beam 
efficiency of 0.74. Data from the WBS spectrometer are used 
and the H and V polarizations have been averaged after visual 
inspection. Table [1] gives the noise level of the observation T^^s 
for the nominal resolution of the WBS spectrometer (1.1 MHz). 
Data were reduced with the Herschel interactive processing 
environment (HIPE) package version 2.9. 



molecular column density of the upper or lower level A^"pp'=' and 

flower ^y 



3. Results 

Spectra of detected target lines are given in Fig.[T] those without 
detection in the Appendix lAl We have detected the ground state 
lines of CH, CH+, OH+, H2O+ and NH. The Unes of OH+, H2O+ 
and NH are seen in absorption without clear signs of emission on 
top of the absorption. The CH line is almost purely in emission. 
The CH^(7 =1-0) line is in absorption, with emission at the 
systemic velocity of the envelope. The CH^(y = 2-1) line is in 
emission. Except for CH^(y = 2 - 1), no line between excited 
states is detected. Not detected are SH^ and NH\ While SH+ 
has been seen in other regions dBenz et al]|2010t iMenten et alj 
[2010), interstellar NH^ has not yet been found. 

By their velocity shift, the absorption lines of CH, CH^, NH, 
OH^ and HaO^ can be roughly classified as (i) a component 
around visr ~ km s ', red-shifted with respe ct to the systemic 
velocity of the envelope at Vsys = -5.5 km s"^ dvan der Tak et alj 

|1999|) and (ii) a blue-shifted component at 11 km s~^ (NH) 

or ~ -16 km s"' (OH+, H2O+ and CH+). 

The velocity integrated emission or absorption is given in 
Table [Tj The absorption per frequency is r = -log(riine/r(;ont), 
with the continuum intensity Tcont in a single side band. The 
velocity integrated emission and absorption are related to the 
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with the Einstein-A coefficient A„; and the frequency v„/ of 
the transition and the statistical weights gu and gi of the lev- 
els. These equations assume optically thin emission or absorp- 
tion without re-emission. For Table [T] contributions of different 
fine/hyperfine components are summed up, weighted by their 
statistical weights. Detected lines together with the upper limits 
from higher levels yield upper limits on the excitation temper- 
ature and the column density, assuming the excitation tempera- 
ture of all lines to be the same. For CH, we obtain Tex < 22.6 
K and N(CH) < 1.6 x 10'^ cm^, OH+ yields Tex < 13.2 K and 
A^(OH+) < 1.2 X 10'4 cm-2 and H2O+ gives Tex < 18.5 K and 
A^(H20^) < 2.3 X 10'^ cm"^. We refrain from deriving column 
densities for CH^ using this method, as the molecule shows clear 
signs of both emission and absorption. 

3.1. Modeling the line spectra 

The line profile is modeled using a two-layer model (red and 
blue) solving the radiative transfer equation along the line of 
sight (Appendix [B). Slab geometry is assumed and the absorb- 
ing/emitting foreground layers cover the whole background con- 
tinuum. Free parameters per component are the column density, 
the excitation temperature and the line width and position. The 
intrinsic line profile is assumed to be a Gaussian. The upper lim- 
its of the OH^ and H20^ excitation temperature are lower than 
the energy of the upper level and for such low temperatures, the 
absorption lines do not depend much on the excitation temper- 
ature. We thus fix Tex at 3 K. Parameters of the best-fit model 
are given in Table|2| The model spectra are overlayed on the ob- 
served spectra in Figure[T] The figure also gives only one compo- 
nent, showing the necessity of a second component. A good fit is 
obtained for the CH, NH and H2O+ fines. For CH and OH+, fur- 
ther components would improve the fit. For OH^, an asymmetric 
line would presumably be needed to fit observations. However, 
for the conclusions drawn here, we concentrate on the strongest 
components. A third component has only been added to model 
the emission peak of CH^ at the systemic velocity of the en- 
velope. The CH column densities derived with this method are 
larger compared to the calculation assuming optically thin lines. 
This is because the fitting of the different fine structure compo- 
nents requires a peak optical depth of r ~ 2. 

4. Discussion 

The blue-shifted component of OH^, H20^ and CH^ sug- 
gests ffie association with the outflow lobe directed towards us. 
iMitchell et aLJ (IT98"9h have detected 4.7 //m '3CO(v = 0-1) 
absorption lines at visr ~ -28 km s '. They derive a '^CO 
column density of 1.1 x 10'^ cm ~^ and a tempera t ure of 200 
K. With better spectral resolution. Ivan der Tak et all ('1999') find 

Visr 21.5 km s ', in better agreement with our observations. 

Their lines are broad and asymmetric, with strong wings to the 
blue, likely originating in winds. It is thus difficult to determine 
the line width and position. However, the '^CO lines reach the 
continuum at about -45 km s"' and are clearly broader than the 
0H+, H2O+ and CH+ lines. Except for CH+, where the two- 
component model does not fit this component well, there is no 
evidence for asymmetric lines in this blue-shifted component. 
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Fig. 1. HIFI-Spectra of hydrides observed towards AFGL 2591. The green line represents the continuum level used for the cal- 
culation of the velocity integrated emission or absorption. The vertical blue line indicates the systemic velocity of -5.5 km s"' of 
AFGL 2591. The red line shows a fitted model spectrum discussed in Section [TTl with parameters given in Table |2] The blue line 
indicates a model spectrum with only one of the components, showing the relative intensity of the different hyperfine components. 
The spectra are corrected for the beam efficiency and the continuum level is given for one sideband (SSB). 



Table 1. Observed Hnes in AFGL 2591. 



Species/Transition/Electronic state/Line 


components'^ 


Frequency 


E„ 


A,,/ 


n,. 


J T'i'^ 


J '"''v 


Tmis 


Af"Ppei(X) 


^'""'(X) 










[GHzl 


[K] 


[s-'] 


[min] 


[Kkms-'] 


[kms-'] 


[mK] 


[cm-2] 


[cm-2] 


CH 


yf,P = 3/22._-l/2i.+ 


^n 


3 


536.7611 


25.76 


6.4(-4) 


4.2 


3.46±0.04° 




16 


3.0(12) 




CH 


yf,/. = 5/23,+ -3/2,,- 


-n 


3 


1661.1074 


105.48 


3.8(-2) 


20.2 


<0.93 


< 0.34*-' 


220 


< 1.3(11) 


< 2.6(12) 


CH+ 


7=1-0 


'j:+ 


1 


835.1375 


40.08 


6.4(-2) 


40.4 


0.91 ±0.03'' 


32.63 ± 0.06'' 


16 


1.9(11)'' 


9.3(13)'' 


CH+ 


y = 2-i 


'r 


1 


1669.2813 


120.19 


6.1 (-2) 


40.4 


3.7 ±0.2 




115 


3.3(11) 




0-H,0+ 


^KiiK^J = ln.3/2 - O00.1/2 


-Bl 


5 


1115.2041 


53.51 


3.1(-2) 


39.2 




7.65 ± 0.05 


33 




1.6(13) 


0-H2O+ 


NkoK^J = 2o2.3/2 - 1 11,3/2 


-Bl 


7 


746.26 


89.33 


5.5(-4) 


40.6 


<0.16 


< 0.52 


25 


< 2.6(11) 


< 3.0(13) 


0H+ 


Nj,F = 1 1,3/2 -0l,3/2 


'I- 


4 


1033.1186 


49.58 


1.8(-2) 


40.6 




21.98 ±0.06 


32 




1.1(14) 


0H+ 


Nj,F = 2i,3/2 - 1 1,3/2 


's- 


4 


1892.2271 


140.39 


5.9(-2) 


21.0 


<0.83 


< 0.28 


210 


< 8.2(10) 


< 2.5(12) 


NH 


^J.Fx.F - l2,5/2,7/2-0l, 3/2,5/2 


-'I- 


21 


974.4784 


46.77 


6.9(-3) 


7.0 




3.90 ± 0.09 


43 




2.9(13) 


NH 


Nj.Ff.F = ll,3/2,5/2-0l,3/2.5/2 


'I- 


21 


999.9734 


47.99 


5.2(-3) 


38.8 




3.99 ± 0.05 


29 




5.0(13) 


NH+ 


y^ = 3/2_-l/2+ 


^n 


5 


1012.5400 


48.59 


5.4(-2) 


38.8 


<0.16 


<0.19 


30 


< 5.9(9) 


< 1.7(10) 


NH+ 


Jp = 3/2+ - l/2_ 


2n 


7 


1019.2107 


48.91 


5.5(-2) 


40.6 


<0.12 


<0.13 


22 


< 4.4(9) 


< 1.2(10) 


SH+ 


!^J.F = l2,5/2 - 0l,3/2 


'i- 


3 


526.0479 


25.25 


9.6(-4) 


4.2 


< 0.083 


< 0.53 


11 


< 4.7(10) 


< 4.5(12) 



The frequency is given for the strongest fine/hyperfine component. The velocity integrated intensity or absorption is summed over all 
fine/'hyperfine components. A(B) = A x 10^. "Icr = ^hVSvT^s, with AV ~ 10 km s 
(corresponding to 1.1 MHz). ''Icr = -^[W6\'T^JT^„^y 
integration time (on+off). 'Hyperfine or fine structure 



(expected line width) and the channel width 6v 
^Non-detections are given by the 3cr value. ''Emitting or absorbing part only. ''Total 



The component is not clearly detected in CH. The NH absorp- 
tion is less shifted to the blue and narrower. The velocity shift of 

10 km s ' corresponds to a narro w absorption feature see n 

in i2cO(7 = 3-2) JCMT observations (va n der Tak erai1[T999l) . 
They conclude that the absorption must take place in cold and/or 
tenuous gas, as it is not found in '^CO(7 = 6-5). 



The red-shifted absorption at Vki- = -0.2 to 2.8 km s is 
likely associat ed to a cold and extended foreground cloud, ob - 
served in CO (JMitchell et al.lll992l:lHasegawa & Mitchelll995h . 
The small line width in their spectra approximately agrees with 
the naiTow absorption lines seen in CH, OH^, HaO^ and NH. 
The broad absorption of CH^ is surprising: since it is only seen 
in CH+(y = 1-0), but not CH+(y = 2-1), the excitation 
temperature must be lower than 8 K, used to fit the blue shifted 
component. This indicates either less dense or colder gas. The 



ionized species OH^, HiO^ and CH^ are more red-shifted and 
have broader lines, compared to the neutral CH and NH. 

The abundances of molecules in the absorption features are 
roughly estimated from the total column density (A^h = NQri) + 
2A^(H2)) obtained from '^CO column densities and presented in 
Tabled We use a CO abundance of «(C O)/«(H7,) = 3.7 x 10' ^ 
(iDotv et al.ll20"02h and '^CO/'^CO ~ 77 (I Wilson & Roodll 1994 . 
Since the red component is a low-density foreground cloud, 
the fractional abundances of CH, HiO^ and NH have been 
compared to a che mical model for a diffuse/translucent cloud 
(iLe Petit et al.ll2004 ) and found to agree within a factor of ~ 2. 
The underprediction of CH^ by a factor of more than 100 is 
well known for such models, but models of turbulent dissipation 
regions (Godardet al. 2009) can explain fractional abundances 
reached here. OH^ is underpredicted by a factor of 6 but the ob- 
served 0H^/H20^ ratio of 7.3 is comparable to the value of > 4 
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found by iGerin et aP (l2010l) and explained by low density gas 
with a low molecular fraction. 

The blue component is thought to arise in tenuous warm 
gas in the outflow lobes directed toward us. The observed OH^ 
and H20^ abundances are lower by factors of 7 and 16, respec- 
tively, but the 0H"^/H20^ ratio i s still large, about 3. As shown 
in Figure 2 of lGerin et alJ (1201 Ol) such ratios are consistent with 
models with densities up to a few thousand cm"^ and enhanced 
radiation fields, as expected in the outflow lobe. The fractional 
abundances of CH^ and NH are similar to those found in the 
red component within a factor of ~ 2. More specific modeling 
of these components must await better determinations of their 
physical conditions. 

Where does the CH and CH^ emission emerge from? The 
line velocity suggests that the emission stems from the envelope 
of AFGL 2591. Abundances of CH and CH+ are difficult to es- 
timate from column densities, due to the strong gradient in den- 
sity in the envelope and thus the changing ex citation. A spheri- 
cal po wer-law density profile of the envelope (Ivan der Tak et alJ 
1 1999 1) yields a density of n(H2) ~ 10* cm -^ at scales of the 
Herschel beam. Relative to the column density of A^h ~ lO^"* 
cm"^ within the region with density > 10* cm"-', we get abun- 
dances of 4 X lO'^ (CH) and 9 x 10" (CH+). 

Why are the hydrides mostly in the ground-state? The exci- 
tation temperature Tgx of a species, assumed to have only two 
levels, is given by 



Table 2. Fitted column density, excitation temperature, 
width and position per component (see Section lTTl i. 



line 



T^kir 



(kTUAE) ■ ln(l +^ «„i,/«(H2)) + 1 



(2) 



with the kinetic temperature of the collision partner Tkin, the en- 
ergy difference between the levels AE, the escape probability 
for a line photon yS, the critical density «ciit and the density of 
a collision partner «(H2). The observed excitation temperature 
of CH+(y = 1 - 0) of 38.4 K can be reproduced by Eq. |2] with 
jS ~ 0.04, assuming a kinetic temperature of 500 K, necessary 
for the formation of the molecule, and a H2 density of order 
10* cm"-', typical at scales of the Herschel beam. This /3 corre- 
sponds to an optical depth of a few at the line center, depending 
on the expression used for /3 = /3{t). Far infrared radiation by 
the dust continuum also pumps the molecule. Critical densities 
of other hydrides are not well known, but likely of similar order. 
Thus, the observed low excitation temperature found for OH^ 
and H20^ can be well explained by the origin of the absorption 
in tenuous gas. 

Predictions for CH and CH^ abundances and line fluxes 
can be made utilizing the photochemical and outflow model 
of [B ruderer et al. (2009b). In a spherical envelope model, 
iBruderer et al.l (12010) find volume averaged abundances orders 
of magnitude lower than observed (~ 6 x 10"" for CH and 
~ 5 X 10"'* for CH^). They present also a model including a cav- 
ity etched out by the outflow, where protostellar FUV radiation 
can escape and directly irradiate the outflow wall. The central 
protostar with a luminosity of Lboi ~ 2x 10"* L© and T^ff ~ 3x10'* 
K can heat the gas in the outflow walls to above 1000 K at scales 
of the Herschel beam. Depending on the radius over which the 
abundance is averaged, abundances of (0.6 - 2) x 10"*^ for CH 
and (0.4 - 4) x 10"^ for CH+ are found with the outflow wall 
model, in acceptable agreement with the observations. 

Calculating the radiative transfer in the fines for direct com- 
parison with observations yields fluxes of (4.3 - 16.1) K km s"' 
forCH+(y = l-0)and(8.9-30.8)Kkms"i forCH+(y = 2-1). 
Note that this calculation includes the influence of the pumping 
by far IR radiation. In the view of uncertainties (such as distance, 
inclination, Tim at formation, and additional contribution to CH^ 



Species 


Component 


N(X) 
[cm-2] 


[K] 


Av 
[kms-'] 


Vkr 

[kms-'] 


A'(X)/A'h 


CH 
CH+ 


CENTER 


4.0(14) 
8.5(12) 


6.6 
38.4 


4.4 
7.0 


-5.9 
-4.6 


4(-9)'' 
9(-ll)'' 



CH 
CH+ 

H2O+ 

0H+ 

NH 974 GHz 

NH 999 GHz 



RED 



2.6(14) 
1.2(14) 
8.3(12) 
6.1(13) 
1.5(13) 
1.5(13) 



1.8 

15.0 

4.5 

5.1 

1.5 

1.5 



-0.2 

1.6 

2.2 

2.8 

-0.2 

-0.1 



2(-8)' 
l(-8)' 
7(-10)^ 
5(-9)' 
l(-9)' 
\(-9y 



CH+ 


BLUE 


1.8(14) 


8.0 


6.0 


-17.1 


4(-9)'' 


H,0+ 




5.4(12) 




13.5 


-16.4 


I(-IO)'' 


0H+ 




1.6(13) 




12.4 


-16.1 


3 (-10)'' 


NH 974 GHz 




1.9(13) 




5.7 


-12.4 


4(-10)'' 


NH 999 GHz 




3.8(13) 




4.0 


-10.4 


8(-10)'' 



Fractional abundances are given relative to A'h = A'(H) + 2A'(H2). 
A{B) = Ax 10^. "fixed at 3 K. '' using Nn ~ 10^^ c m"^ see text. 
' using A'h ~ 1.2 x lO^^ cm"^ ( iMinh & YanlpOSh . 
'' using A'h ~ 4.6 x 10^^ cm"^ ( JMitchell et alJl98gh . 



by rea ctions of C^ with vibrationally excited H^i lAgundez et al.) 
2010j, we consider the agreement with observations good. The 
abundances of OH+ and H20^ are also enhanced in the outflow 
walls, but lack of collisional rate coefficients prevents us from 
making quantitative model results. 

We conclude that the scenario of FUV irradiated outflow 
walls can quantitatively reproduce the abundances and fluxes ob- 
served in CH and CH^ emission. This indicates the presence of 
a large volume/mass of hot, FUV irradiated and extended gas in 
the envelope of this high-mass young stellar object. 

5. Conclusion 

We use the new Herschel-HWl instrument to observe light hy- 
drides towards the high-mass star forming region AFGL 2591. 
Besides OH^, H20^ and NH detected in absorption only, we see 
CH and CH^ in emission. The red and blue velocity shifted ab- 
sorption can be assigned to a foreground cloud and an outflow 
lobe directed towards us, respectively. The emission of CH and 
CH^ is at the source velocity. Its abundance and excitation can 
be explained in the scenario of direct FUV irradiation of the en- 
velope through a low density cavity of the outflow region. 
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Appendix A: Spectra of non-detections 



Figure lATTI presents spectra without detected target lines, given 
in Table [U For clarity, the spectra of CH(5/23,+ - 3/22,-) and 
OH^(2i3/2 - 1 1,3/2) have been scaled by a factor of 0.1. In the 
spectrum of CH(5/23,+ - 3/22,-), the H20(22i - 2i2) line is 
detected, which will be presented in an upcoming paper (van 
der Tak et al., in preparation). The H20^(2o2,3/2 - 1 11,3/2) spec- 
trum contains a S O2 line with uppe r level energy similar to the 
lines detected by van der Tak et al.l (12003) in the same source. 
In the other sideband of the NH'^(3/2+ - 1/2-) spectrum, the 
OH^(2i,3/2 - 1 1,3/2) line is seen. In the upper side band of the 
SH+(l2,5/2 - 0i,3/2) hne, the HCO+(6 - 5) line has been detected. 
The velocity range of ±40 km s ' corresponds to a frequency 
range of ±67 MHz at 500 GHz and ±253 MHz at 1900 GHz. 
Thus, the hnes should be seen in the spectra, even if the pre- 
dicted rest frequency is grossly wrong. 
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NH + (3/2^-1/2_) 

HCD^(6-5) in USB 
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^- ^^^ I ^^^^^^^^ ■^^^"^A.^^ ^a 
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fine/fine components. The line opacity is calculated from the 
sum over the M fine/hyperfine components of the transition. 



r'iv) = A^; 



mol 



z 



" 4* 

87r(v*)2 "' 



81 



¥Hv) . 



(B.2) 



with the column density per layer N'^^^ the Einstein- A coeffi- 
cient A^^p the normalized line profile function </i''*(v), the normal- 
ized level population of the upper and lower level (x|; and x'^ , 
respectively) and the statistical weights {gu and gi, respectively). 
The normalized level population per layer is a Boltzmann distri- 
bution for the temperature T'^^. The line profile function <p''''(v) 
is assumed to be a Gaussian with width Av' (FWHM) centered 
at the frequency of the component minus the velocity vj^.^. of the 
layer. 

The free parameters per layer are thus the excitation temper- 
ature Tg^, the column density A^^^], the width Av' and the veloc- 
ity vj^j.. Two layers are used for CH, H20^, OH^ and NH with 
a red-shifted layer in front of a blue-shifted layer. Three layers 
are used for CH^ with a layer centered at the source velocity 
between the background continuum and the blue-shifted layer. 
To constrain the parameters of the layer, the least squares of the 
modeled to the observed spectra are minimized. 

Molecular data of the individual hyperfine components of 
CH, OH^, NH and H20^ used to model the line spectra in Figure 
[T]are given in Table IB. II The table gives the frequency, Einstein- 
A coefficient (A,,/), the upper level energy (£„) and the statistical 
weights of the upper and lower level of the transition (gu and 
gi, respectively). The third column refers to the velocity shift 
in km s ' of the component relative to the one with the largest 
Einstein-A coefficient. The component with the largest Einstein- 
A coefficient is marked by (*). 



Fig.A.l. Spectra observed towards AFGL 2591 without detec- 
tion of target lines. The velocity is given relative to the systemic 
velocity Vsys = -5.5 km s"'. 



Appendix B: IVIodeling the line spectra 

This appendix gives the equations used to model the line spec- 
tra. The intensity t'"°'^''^{v) [K] at frequency v for a slab of A^ 
emitting or absorbing layers in front of a continuum source is 
obtained from 



(V) = 7^contmuum(^y-Ef,,AT'(y) ^g j^ 



T^model^,,\ _ ^continuum 

N 2 

C 



,.1 2v?,^ 



with the (background) continuum intensity r'^°"'(v), the opac- 
ity per layer At'(v), the frequency of the component vo with 
the largest Einstein-A coefficient and the Planck function By(T). 
Note that this equation assumes a covering factor of 100 % for 
each layer and the same excitation temperature for all hyper- 
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Table B.l. Hyperfine components used to model the spectra. 



Species 



Frequency 
[GHz] 



Siiift to (*) 
[km s-'] 



Eu 
[K] 



[s-i] 



Si 



CH(*) 


536.7611 


0.0 


25.76 


6.4(-4) 


5 


3 


CH 


536.7820 


-11.6 


25.76 


2.1 (-4) 


3 


3 


CH 


536.7957 


-19.3 


25.76 


4.3(-4) 


3 


1 


H2O+ 


1115.1505 


14.4 


53.51 


1.7(-2) 


4 


2 


H2O+ 


1115.1862 


4.8 


53.51 


2.8(-2) 


2 


2 


HzO^ (*) 


1115.2041 


0.0 


53.51 


3.1(-2) 


6 


4 


H2O+ 


1115.2632 


-15.9 


53.51 


1.4(-2) 


4 


4 


H2O+ 


1115.2989 


-25.5 


53.51 


3.5(-3) 


2 


4 


OH+ 


1032.9979 


35.0 


49.58 


1.4(-2) 


2 


2 


OH+ 


1033.0044 


33.1 


49.58 


3.5(-3) 


4 


2 


OH+ 


1033.1118 


2.0 


49.58 


7.0(-3) 


2 


4 


OH+ (*) 


1033.1186 


0.0 


49.58 


1.8(-2) 


4 


4 


NH 


974.3156 


50.1 


46.77 


1.8(-6) 


4 


4 


NH 


974.3426 


41.8 


46.77 


6.0(-6) 


6 


4 


NH 


974.3546 


38.1 


46.77 


4.5(-6) 


4 


2 


NH 


974.4106 


20.9 


46.77 


8.9(-5) 


4 


6 


NH 


974.4114 


20.6 


46.77 


5.1 (-4) 


2 


4 


NH 


974.4363 


12.9 


46.77 


2.3(-3) 


4 


4 


NH 


974.4375 


12.6 


46.77 


1.3(-3) 


6 


6 


NH 


974.4440 


10.6 


46.77 


1.8(-3) 


4 


4 


NH 


974.4504 


8.6 


46.77 


4.8(-3) 


2 


2 


NH 


974.4622 


5.0 


46.77 


5.0(-3) 


4 


2 


NH 


974.4710 


2.3 


46.77 


5.7(-3) 


6 


4 


NH 


974.4754 


0.9 


46.77 


3.4(-3) 


4 


2 


NH(*) 


974.4784 


0.0 


46.77 


6.9(-3) 


8 


6 


NH 


974.4793 


-0.3 


46.77 


6.0(-3) 


6 


4 


NH 


974.5313 


-16.3 


46.77 


2.6(-4) 


4 


6 


NH 


974.5398 


-18.9 


46.77 


6.5(-4) 


2 


4 


NH 


974.5581 


-24.5 


46.77 


9.8(-4) 


2 


2 


NH 


974.5648 


-26.6 


46.77 


7.7(-4) 


4 


4 


NH 


974.5744 


-29.5 


46.77 


8.1 (-4) 


6 


6 


NH 


974.5830 


-32.2 


46.77 


2.6(-4) 


4 


2 


NH 


974.6078 


-39.8 


46.77 


1.2(-4) 


6 


4 


NH 


999.8784 


28.5 


47.99 


9.9(-4) 


6 


4 


NH 


999.8821 


27.4 


47.99 


3.9(-3) 


4 


4 


NH 


999.9060 


20.2 


47.99 


3.1(-7) 


4 


4 


NH 


999.9079 


19.6 


47.99 


3.0(-3) 


2 


4 


NH 


999.9192 


16.3 


48.00 


l.l(-3) 


2 


4 


NH 


999.9212 


15.7 


47.99 


9.6(-4) 


4 


2 


NH 


999.9451 


8.5 


47.99 


1.8(-3) 


4 


2 


NH 


999.9470 


7.9 


47.99 


1.8(-4) 


2 


2 


NH 


999.9582 


4.5 


48.00 


1.7(-3) 


2 


2 


NH(*) 


999.9734 


0.0 


47.99 


5.2(-3) 


6 


6 


NH 


999.9771 


-1.1 


47.99 


3.9(-4) 


4 


6 


NH 


1000.0010 


-8.3 


47.99 


2.9(-3) 


4 


6 


NH 


1000.0068 


-10.0 


47.99 


1.2(-3) 


6 


4 


NH 


1000.0105 


-11.1 


47.99 


1.7(-3) 


4 


4 


NH 


1000.0288 


-16.6 


47.99 


4.9(-4) 


4 


2 


NH 


1000.0345 


-18.3 


47.99 


1.6(-3) 


4 


4 


NH 


1000.0363 


-18.9 


47.99 


4.4(-5) 


2 


4 


NH 


1000.0476 


-22.2 


48.00 


4.5(-3) 


2 


4 


NH 


1000.0527 


-23.8 


47.99 


l.l(-3) 


4 


2 


NH 


1000.0546 


-24.3 


47.99 


4.1(-3) 


2 


2 


NH 


1000.0659 


-27.7 


48.00 


1.6(-4) 


2 


2 



A(B) = A X 10'^. 



